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Abstract 



We study the Infra Red (IR) limit of dispersion relations for scalar and 
fermion fields in a Lifshitz-type Yukawa model, after dressing by quantum 
fluctuations. Relativistic-like dispersion relations emerge dynamically in the 
IR regime of the model, after quantum corrections are taken into account. 
In this regime. Dynamical mass generation also takes place, but in such a 
way that the particle excitations remain massive, even if the bare masses 
vanish. The group velocities of the corresponding massive particles of course 
are smaller than the speed of light, in a way consistent with the IR regime 
where the analysis is performed. We also comment on possible extensions of 
the model where the fermions are coupled to an Abelian gauge field. 



1 Introduction 

Quantum field theories in the Lifshitz context exhibit interesting renormaliza- 
tion properties pQ. Lifshitz type models are based on an anisotropy between 
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space and time directions, which is characterized by the dynamical critical 
exponent z, determining the properties of space-time coordinates under scale 
transformations t — )■ b^t and x — )■ fex . For z > 1 the higher powers of momen- 
tum in the propagators lower the superficial degree of divergence of graphs, 
yielding the renormalizability of new interactions, such as the four-fermion 
interaction [2j. In the specific case of z = 3 in 3+1 dimensions, the scalar 
field is dimensionless, so that any function of the field is renormalizable. 
The case of an exponential potential for the scalar field has been studied in 

where the exact effective potential has been derived and found consis- 
tent with the well-known Liouville theory in 1+1 dimensions. Moreover, in 
this framework, divergences of renormalizable interactions in the Standard 
Model become softer, for instance in the Yukawa model [1] only logarithmic 
divergences appear. Finally, we mention that this approach has also been 
proposed as a renormalizable alternative to Einstein General Relativity [5], 
which could lead to a consistent model of Quantum Gravity. 

While absent from the classical action, Lorentz symmetry can naturally 
be generated in Lifshitz-type models through quantum corrections, since the 
corresponding kinetic term is a relevant operator and dominates the disper- 
sion relation of the modes in the infrared regime (IR). In [6j, the authors 
study the case z = 2, for scalar theories involving derivative interactions. 
There, logarithmic divergences imply a speed of light running with a mass 
scale, and the discrepancy in the IR speeds of light which are obtained for 
different species of particles is discussed. 

In the current work we shall study a softer version of the Lifshitz frame- 
work, namely a. z = 3 Lifshitz-type Yukawa model, where the only divergence 
appears in the scalar mass. As a consequence, there is no running speed of 
light in the model, while the pertinent IR dispersion relations are consistent 
between different species of particles. However, an important feature of the 
model, which will be discussed in section 2, is the fact that the quadratic 
term in the scalar- field derivatives has the "wrong sign" , so that an appar- 
ent tachyonic mode develops in this theory. Nevertheless, as a consequence 
of the six-order derivative term, which has the correct sign, the tachyonic 
mode disappears and a relativistic dispersion relation can be obtained after 
expanding the frequency around its minimum. This is similar to an expan- 
sion around the Fermi surface in condensed matter physics, and leads to a 
momentum shift. 

The structure of the article is as follows: in section 2, we discuss the fate 
of a tachyonic mode to sixth order in derivatives of the scalar field, and the 
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recovery of a relativistic-like dispersion relation for the relevant excitations 
in the model. The corresponding effective light cone depends on the particle 
species, but we explain that no contradiction with the speed of light arises. 
In section 3, we study the non-perturbatively dressed masses of the system, 
and show that these do not vanish, even if the bare masses go to zero. We 
show this by means of a non-perturbative treatment of the model, in the 
form of a differential Schwinger-Dyson approach, based on a set of self con- 
sistent coupled equations for the dressed parameters. We exhibit numerical 
solutions for these self consistent equations, where one can see that the no 
IR divergence occurs, unlike the case of the perturbative analysis based on a 
massless bare theory. 

In section 4, we discuss an extension of the model where the fermions 
couple to an Abelian Gauge field (a form of anisotropic Quantum Electro- 
dynamics - QED). We show that the one-loop analysis is not sufficient to 
allow for an interpretation of the gauge boson as a photon in such a model, 
in the sense that its group velocity is different from that of light in vacuo. 
Conclusions and Outlook are presented in section 5, where we stress that the 
dynamical mass generation in our Yukawa-Lifshitz model allows fermions to 
have a mass, without Higgs mechanism, and that the possibility of formu- 
lating a Lifshitz Standard Model may lead to phenomenologically relevant 
realizations. 

2 IR dispersion relations 

We explain in this section the different features related to the IR of the 
dressed theory, starting from a Lifshitz-type Yukawa model. For the sake of 
clarity, the corresponding calculations can be found in Appendix A. 



The Lifshitz-type Yukawa model in D + 1 dimensions and with anisotropic 



2.1 The Model 



scaling z is defined by 





(1) 



Throughout this work we use units where the speed of hght in vacuo is c = 1. 
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where A = d^d^Sij, with i,j spatial indices. We consider here the case 
z = D = 3, where the mass dimensions of coordinates and fields are 

[t] = -3 [x] = -l [<P]=0 [^] = | (2) 
The bare action is then 



2'"' 2 

^o<P'^ - mlipip - gocpipip I , (3) 



where the Yukawa coupling qq has mass dimension [qq] = 3 and mo is the 
bare mass of particles. We consider the same bare mass for the scalar and 
the fermions: a different mass would involve a new parameter which would 
not change qualitatively the non-perturbative analysis we present in the next 
section. This theory is super-renormalizable, and the only divergence occurs 
in the dressed scalar mass 



2.2 Naive approach 

Quantum fiuctuations generate lower orders in space derivatives of the field, 
such that the effective theory is 

where M and m are respectively the dressed scalar and fermion masses, 
and X^, A^, C,/,, Cv^ arise from quantum corrections. Sint contains interactions, 
higher orders derivative terms, derivative interactions and higher orders in 
fermion fields. 

If one considers the free bosonic and fermionic sectors individually, then there 
are several possible ways to rescale space time coordinates and fields, in order 
to recover the relativistic dispersion relations. But because of the interac- 
tions contained in Smt, the space time coordinates have to be rescaled in 
the same way for the bosonic and the fermionic sectors. Indeed, the Yukawa 
interaction, as well as other interactions generated by quantum fiuctuations. 
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are local and occur at the same event in space time. As a consequence, we 
define the global rescaling 

t ^ at X -> 6x (j) ^ Acj) ip Bip, (5) 

where the mass dimensions are [b] = [B] = 0, [a] = —2 and [A] = — 1, in such 
a way that the new mass dimensions of space-time coordinates and fields are 
the expected ones in an isotropic theory. The new quadratic terms of the 
effective action, then, become: 

ab' I dtd'^i^A^ (^('^)' - ■ (9<P) - ^0^) (6) 

\ a b 

It is easy to check that one cannot set all the coefficients of the derivative 
terms to unity in a consistent way. Nevertheless, one can impose the require- 
ment that the coefficients of the time derivatives be equal to unity, which 
leads to the following constraints: 

(1 + C^)^'^' = a and (1 + C^)fi'6' = 1 . (7) 

We find that no further constraint can be imposed in a consistent way. We 
are then led to the following quadratic terms in the effective action 

+ i^ljdo'y^^l) — iv^'il){'y ■ d)-ip — m%l)'ijj^, (8) 

where we defined 

= X^A^ab , = X^B'^ab'^ 

M2 = M^A^ab^ , m = m^B^ab\ (9) 

Note that the constraints ([7j) do not fix the parameters a, b, A, i? in a unique 
way. The new parameters (Q (which are not independent) are renormalized 
parameters, fixed in principle by the "experiment", and define a specific 
rescaling, with unique values for a, b, A, B. 
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2.3 Non-trivial scalar dispersion relation 

The previous arguments hold if A,^ > 0. However, as demonstrated in Ap- 
pendix A, the one-loop expressions for X^f, and in the model are 

Ai^^ = ^ < and A« = -|4^ (10) 



with 



dxx' /8 3(3x^ + 7) 15x' \ 

' (1+X6)3/2V3 4(l + x6) 2(l + x6)2y' 

J = I ^^^^0.16 (11) 



/o (1 + a;6)5/2 

As a consequence, the scalar dispersion relation is of the form (before rescal- 
ing fields and space time coordinates) 

(1 + Qco^ = M'- \X^\e + r^^k"" + (1 + ^^)k' + 0{k'), (12) 

where ?]<^ and are quantum corrections, such that 

ri^ocgl/ml and ^^(xgl/ml. (13) 

This dispersion relation presents an apparent tachyonic sign, but because of 
the positive term fc®, the frequency has actually a non-trivial minimum for 
some momentum ko such that 



iaW' 



kt = \^,-\^, = \l'-^ + 0{gt/ml). (14) 



One can therefore expand up in powers of q = k — k( 

o o 1 (fu)"^ 



q\ (15) 

ko 



2 dk 

and we obtain from eq.( !T^ 

uj^ ~ /i6 + 4|Aj)|g2 + 0(Vmo)' (16) 



with /i« = M** - 3 I = + C(^o/mo) 



^One can recall at this stage a similar situation that characterizes condensed- matter 
systems, linearized about their Fermi surfaces. 
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such that this momentum shift imphes the replacement of A^^^ by 4|A^^''| in 
the new dispersion relation, based on the shifted momentum q. We check 
now that this shifted momentum cannot be detected at one loop (order Qq). 
Indeed, the corresponding kinetic term in the action is, in Fourier compo- 
nents, 

ko)^ - /^^)0<^,k0-a; -k (17) 
j 0w,k+ko</'-w,-k-ko 

M6)0^,u0-^,_u + C((7o/mo)', 

since A</, = 0{gQ/ml) and /cq = O^go/mo), and the linear term in ko is a 
surface term. From eq. (|T71) . the scalar dispersion relation is, at one loop, 

u;2~M^ + 4|a5,'^|A;2 (18) 

We finally note that the fermionic dispersion relation always involves A^, and 
is therefore independent of the sign of A^. 



2 J 271 (27r)3 V ' 



2.4 Group velocities and effective Yukawa coupling 

From the quadratic action ([8]), we can obtain the IR dispersion relations for 
the fermionic or bosonic excitations in the model 

U;2 = ^2 ^ ^2^2 ^ 0{kl^l)\ (19) 

where f is a generic speed {v^ or w^) and yU is a generic dressed mass (m 
or M). Note that the speed v cannot be identified with the speed of light. 
Indeed, the speed of light would be obtained from the dispersion relation 
( !T9|) in the limit where /i — )■ 0, but we will show in the next section that the 
scalar and fermion dressed masses never vanish, even if the bare mass mo 
goes to zero, as a consequence of dynamical mass generation. Instead, one 
can find from the relation f lT9|) the product of the group velocity Vg and the 
phase velocity Vp of particles: 

vgv, = ^r = -" + ^(k/f^r . (20) 
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A word of caution is in order at this point. One should keep in mind that 
this dispersion relation is valid for k << fx only. In this regime of momenta, 
V cannot represent a limiting speed of particles in the model, since the latter 
can only be obtained in the limit k ^ oo. 

From the constraints ([7]) and the definitions (Q , we find the following relation 

between and 

Using the one-loop expressions ffTOl) . we can readily see that the speeds 
and V(f) satisfy the relation 

^ ~ ^ ^ ^ « 1 (22) 

2 /lA^^^I 327rv/7 "^o 



Finally, one can express the dressed dimensionless Yukawa coupling g in 
terms of the different parameters obtained after the rescaling ([5]). From the 
interaction 

g H^^jj = g ab^AB^ (fi^i/j, (23) 
and the relations (J7f9i) . we obtain 

We show in Appendix A that the one loop coupling g^^^ is equal to the bare 
coupling go, so that the regime remains perturbative at one loop {g « 1), 
if goV^J"^ « |4A^^^|^/^ or, equivalently: 

rrj 1/3 1/3 

V, « ^ ^ 0.31 ^ (25) 

V vr^ mo mo 

It is interesting to observe f l22|) that the speed of the fermions is smaller 
than the corresponding one, f^, for the scalar fields. However, the reader 
should recall (c.f. fl20|) ) that t>^ and f| represent the respective products 
of the group and the phase velocities and thus should not be confused with 
the propagation velocity of the physical excitations. Nevertheless, in the 
infrared, the corresponding group velocities are: 

vl = -j^= ^ v'y (26) 
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since ^ v'^ . As we shall see in the next session, the effective masses /x^ 

of fermions and scalars are of the same order, and thus, for a given momen- 
tum scale, on account of fj22|) . the group and phase velocities of the fermions 
will be smaller than the corresponding ones of the scalars. Moreover, the in- 
equalities fl22|) and fl25|) show that the model can represent "slowly" moving 
particles only, in the perturbative regime jra^ << 1. 
Finally, the relativistic-like dispersion relation ( I19p exhibits a different effec- 
tive light cone for the scalar and the fermion (although we stress again that 
the speed of light is not questioned by our analysis), showing that this model 
can be relevant to Condensed Matter, as we will discuss in the conclusion. 

3 Non-perturbatively dressed masses 

In this section we shall demonstrate that the dressed masses of the model 
cannot vanish, even if the bare masses do. This dynamical generation can be 
studied only through a non-perturbative approach (such as the Schwinger- 
Dyson approach), since the resulting masses are not analytic functions of the 
coupling constant, and no perturbative expansion can exhibit the dynamical 
mass-generation mechanism. 

We will quantize the theory (|3]), and derive an exact self-consistent equation 
for the proper graphs generator functional F, in the form of a differential 
Schwinger- Dyson equation, giving the evolution of F with the amplitude of 
the bare mass mg. As explained in the original paper [7], this approach con- 
sists in controlling the amplitude of quantum fluctuations: if one starts from 
a large mass mo, quantum fluctuations are "frozen" , and the system almost 
remains classical, F ~ S*. As the bare parameter mo decreases, quantum 
fluctuations gradually dress the system, and it can be shown ^ that the 
corresponding flows in mg are equivalent, at one-loop, to those given by the 
standard Callan-Symanzik equations. 

We first review the path integral quantization of the model, in order to define 
our notations, and then derive the exact evolution equation for the proper 
graph generating functional F with the bare mass mg. Then, upon assuming 
a specific functional form for F, we derive the corresponding dressed scalar 
and fermion masses. The pertinent evolution equations are non-perturbative. 
They constitute a set of self consistent coupled equations, thereby represent- 
ing a resummation of graphs. 
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3.1 Path integral quantization 

The quantum theory is based on the partition function 

Z\j, V,r]] ^ J T^[<f>, i^, V'] exp (^iS + i J dtd'xij^) + rjiP + i^r))^ 

= exp{tW[j,7],r]]), (27) 

where j, 77, r/ are the sources, and W is the connected graphs generator func- 
tional. The classical fields are defined as 

—i6Z —iSZ — i 6Z , . 



and we have 



SjSj 
SrjSf] 



i4'ci4'ci - i <(t)(t)> 

H^d^ci -i<iiip>, (29) 



where 



<...>= 1 y x>[(f), ip, V^](- • •) exp (^iS + i J dtd^x{j(p + rjip + ipr))^ (30) 

The proper graphs generator functional r[0ci, V'd) V'd] defined as the Leg- 
endre transform of W 

r^w- J dtd^x{j(f)ci + v^ci + ^civ) , (31) 

where the sources must be viewed as functionals of the classical fields, and 
therefore functions of the bare parameters in the model. From this definition, 
we obtain 

. sr _ sr 

0(Pcl Olpcl Oljj^ 

and 

(5'r),,,, = -{3'w)jl, {s'r)^,^^^ = -is'w)ri, (33) 

where the notation {S^A)ij represents the (i, j)-element of the matrix with 
components equal to the second derivatives of A. 
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3.2 Exact evolution equation 



From these definitions and properties, we can now derive the exact evolution 
equation of F with mo. In what follows, we denote a derivative with respect 
to mo with a dot. The first step is to note that 

- j dtd^x {dmoJ(pci + fli^d + ^clV) = W (34) 
Using the identities fl28|) to f p3|) . then, we obtain 

f = -3mo J dtd^x (ml < (p"^ > + < tpip >) 

= -3ml j dtd-^x (mlcpli + ip^ipd) 

+3^-oTr {m3(5^F)-,^^ + (5^F)-^ J . (35) 

We stress that the self consistent evolution equation (1351) for F is exact, and 
not based on any assumption or expansion. The resummation of all the 
quantum corrections to the bare action is contained in the trace. Note that, 
if one replaces F by the bare action S in this trace, we obtain the usual 
expression for the one- loop theory, after integration over mo. 



3.3 Ansatz for the proper graph generating functional 

From now on, we omit the subscript cl on the classical fields. In order to get 
physical insight on the quantum theory obtained from the solution of eq.(l35l). 
we need to assume a functional form for F. We consider the following ansatz 

F = y dtd^x!^^{(pf - ^((90) ■ (A^S^) + #7°V' - ^^A(7 ■ d)ij 

-V(0)-f/(0)V^^|, (36) 

where U{(j)),V{(f)) are scalar potentials to be determined. The ansatz fl36!) 
ignores lower order kinetic terms generated by quantum fluctuations, since 
these are negligible compared to in the UV regime relevant for the loop 
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integral (trace in the evolution equation f l35p ). Also, in the framework of 
the gradient expansion, we neglect higher order derivatives, which are also 
generated by quantum corrections. As a consequence, we concentrate on the 
scalar potential sector (self coupling potential V{(j))), and the scalar-fermion 
coupling U{(f))ipip. 

We show in Appendix B that, after plugging the ansatz fl5B]) into the 
exact evolution equation f l35p . the evolution equations for the potentials are 

o 5^2 ml/2A^\ mlU ^ f2A^ 
V = 3m^02 + -|ln ^ +-^ln 



IT \VV" J 7C \ U 



U = 3ml + ^ (^7^^ In (^) {V" + - 2mlU) (37) 



ml mlU" + 2[U'f ^ ml U{UU" + 2[U'' 



8vr2 V'-W^ 87r2 V"{V"-U^) 

where a prime denotes a derivative with respect to the scalar field. Note that 
the apparent singularities in the evolution of U, when U"^ — )■ V", actually 
cancel each other. We indeed checked that, after setting V" = U^{1 + e), we 
obtain 

£> = - ^ + 0(<), (38) 
SO that the evolution equation is well defined. 

3.4 Polynomial expansion of the potentials 

We consider the following expansion for the potentials 

f/(0) = m^ + g^ + 0{^^), (39) 

where m, M are respectively the fermion and scalar dressed masses, and g 
is the dressed Yukawa coupling. The linear term k0 is purely generated by 
quantum fluctuations, and arises from the trace in eq. (l35ll which contains 



(5^r);^^. The field- independent term Vq takes care of the divergences present 
in the evolution equation for V{(f)), since setting = in the latter equation 
gives 

m^, /2A^\ ml . /2A^\ ,^ , 
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and the evolution equations for the remaining parameters g, m, M, obtained 
by differentiation of eqs. fl37p . become then divergence-free. Only the evolu- 
tion equation for k, still contains a divergence: 



7r2 



In 



2A3 



(41) 



However, since this parameter does not enter into the other evolution equa- 
tions for g, m, M, the latter have cut-off independent flows in rriQ. The param- 
eter K decouples from the rest and does not appear in any loop calculation, 
since the linear term in is not present in the bare theory. 
Plugging the expansions fl39|) into the evolutions equations f l37|) . we obtain a 
set of non-linear coupled differential equations for g, m, M: 



5 9'^ ^0 
2 g'^ml + 



(M' 



In 



9' 



rrin 



127r2 M6 - m6 
3g^ml r 



m 
M6 



6\2 



27r2 

9' 



3m^ — rriQ 



M 
m 



(M6 



6^2 



+ 



4m^(m^ — mg) 
(M6 - 



6^3 



47r2 (M6 



m 



6\2 



m'r 



m 
M6 



In 



(42) 



+ 3(mo - m^) 



Once again, one can check that the apparent singularities when — )■ 



actually do not occur: if one sets 



m 



'1 + e), an expansion in e shows 



that the singularities cancel out in the limit e — 0. 

As expected, one can also check that the one-loop results (calculated in 
Appendix A) are recovered after integration over mg, if one replaces the 
dressed parameters in the right-hand side of the evolution equations f H2]) by 
the bare ones: 



Ml, 



m 



(1) 
.(1) 



mo/ 3 



1 



247r2m|] 



90 



(43) 
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where we note that the one-loop correction to the couphng vanishes. 
In figs. 1 and 2 we plot the numerical solutions of eqs. fl42p in the range of 
bare masses from mo/ A = 1 down to mo = 0, for go/ = 0.01. We observe 
that the dressed masses do not vanish when the bare mass mo goes to zero, 
due to dynamical mass generation. This is an important physical feature of 
the model. 

Callan-Symanziknon-perturbativdIowforM Callan-Symanziknon-perturbativ^lowform 




m_0 

Figure 1: Non-perturbative flow for the dressed masses M (scalar), m 
(fermion) versus mo. The dressed masses do not vanish when the bare mass 
mo goes to zero, due to dynamical mass generation. 



3.5 IR stability 

One can see from the one-loop expressions fH3l) that the limit mo — )■ leads to 
IR divergences. This is actually the case at any order of the loop expansion. 
Our numerical analysis shows that the resummation provided by the set of 
the coupled equations fH21) restores the IR stability of the system, and that 
no divergence occurs in the dressed system, when mo — )■ 0. Figs. 3 and 4 
compare the one-loop and non-perturbative flows of M and m in the range 
of bare masses from mo / A = 1 down to mo = (the region around mo = is 
zoomed in, to highlight the difference). We clearly see with the fermion mass 
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Callan-Symanziknon-perturbativdlowforg 
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1.00003E-2 

g 

1 .00002E-2 
1.00001 E-2 

1 E- 2 I — 

0.2 0.4 0.6 0.8 1 

m_0 

Figure 2: Non-perturbative flow for tlie dressed coupling g versus tjiq. Note 
tliat tlie one-loop correction to the coupling vanishes, explaining the small 
change in values for g. 

that the would-be IR divergence, in the limit where itiq — )■ 0, does not occur in 
the non-perturbative flow. This feature is a consequence of dynamical mass 
generation, which was already studied in the context of (2+l)-dimensional 
quantum electrodynamics {QED3), following the present non-perturbative 
method [8]. 

4 Comments on the gauged model 

In this penultimate section, we would like to make some remarks concerning 
a possible extension of this Yukawa model to a gauged one, in which the 
fermions couple to an Abelian Gauge field. First, let us discuss a few features 
of the Lifshitz-type QED, which may be relevant for our purposes. Studies 
with anisotropic higher-order derivatives in QED have already taken place 
in [9], but there, the Lorentz- violating terms are added to the usual Lorentz 
invariant QED. In the current article, we shall discuss a QED-tjpe model, 
where only higher order derivatives are present in the model, as is the case 
of the Lifshitz- Yukawa model studied in previous sections. 
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Callan-Symanziknon-perturbativa/sone-loopflowforM 



Callan-Symanziknon-perturbativa/sone-loopflowforM 




Figure 3: Comparison between the one-loop and the non-perturbative flows 
for the scalar dressed mass M versus mo (the figure on the right is a zoom 
near mo = 0) 



Callan-Symanziknon-perturbativa/sone-loopflowfonn 



Callan-Symanziknon-perturbativa/sone-loopflowform 
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Figure 4: Comparison between the one-loop and the non-pcrturbative flows 
for the fcrmion dressed mass m versus mo (the figure on the right is a zoom 
near mo = 0. 
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We first note that the action for an Abehan gauge field in this framework, 
reads: 

Sa = -\j dtd'^ {2Fo,FO'= + G^.G^'} , (44) 
with k,l = 1, ■ ■ ■ , D and 

Fok = dtA,-d,Ao , Gki = -A{dkAi-diAk), (45) 

where A = 5^^dkdi. The gauge fields have dimensionalities [Aq] = 2 and 
[Ak\ = 0, and the action (jUj) is invariant under the gauge transformation 
A^^A^ + d^e, where [9] = -1. 

In order to find the gauge propagator, we first calculate the second functional 
derivatives of the action (jH]) with respect to the gauge fields. 



- -A 5(t - O'^^'^^x - x') (46) 



6A06A0 



SAkSAi 
In Fourier space, these read 



[r]''\d^ - A') - A^d'd'] Sit - t')5^''\x - x') . 



p2 5(u;-a;')5(^)(p-p') (47) 



5AoSAo 

= cop' Sico ~ u')5^^\p - p') 

= [((py-a;V' + (pWl>T 6{co-u:')6^^\p-p'). 
Hence one can check that, for any fi = 0,1, ■ ■ ■ , D, 

^^Sq 5^Sq S^Sq 

m;m:p^ = u;m,''^sajaJ^ = ' ^''^ 

which shows the absence of an inverse for the operator S'^Sg. 

As usual, this is remedied by the addition of a gauge fixing term. In 
our Lifshitz case, this is naturally provided by the Coulomb gauge, which is 
equivalent to adding the term 

-^A'id'AkY (49) 
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to the Lagrangian, so that the space components of the second derivatives 
(H7|) change to 



= ((p^)=' - u') v"' S{u - uj')6^^\p - p'). (50) 

oAkOAi 

Using the tensorial structure available in the D-dimensional space, we search 
for a gauge field propagator D^^iuj, p) in the form 

Doo = A = Bpk, Dki = Cr]ki + Ep^Ph (51) 

where A,B,C,E are found from the definition 

^D,P = - c-O^^^Hp - P') • (52) 

This leads to the following structures 

^00 = 5fl-7^1 (53) 



Ofe 



p2 V (p2)3 



(P^)^ 



Dkl = T-^v^ Vkl + 



^p2)3 - u;2 ^''^ ■ (p2)4 

In the Yukawa model discussed above, the scalar field is real and has no U{1) 
charge. However, the (Dirac) fermion field can couple to the gauge field. 
This can be done in several ways, respecting gauge invariance. The minimal 
coupling can have one of the following forms 

ip{idt — eAQ)'y'^ip — ip (^id — eA^ ■ ^ ip 

i^{idt - eAo)7V -i^(id- e a) ^ (id - e a) ■ 7 ^ (54) 

Since the gauge coupling is dimensionful ([e] = 1), the quantum theory will 
exhibit fermion dynamical mass generation (the gauge boson remains mass- 
less because of gauge invariance). Therefore the model cannot describe mass- 
less fermions and, as in the Yukawa model, the fermion group velocity will 
be smaller than the speed of light. 

The Lorentz-restoring kinetic terms arising from quantum fiuctuations, for 
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both the fermion and the gauge field, contain several contributions because 
of the new vertices appearing in the minimal couplings flM|) . Nevertheless, 
at one loop, the only contribution to the space components of the polariza- 
tion tensor, which depends on the external momentum, is the same as the 
corresponding one in the Yukawa model, and is proportional to 



As we have already seen in the context of the Yukawa model above, this 
will lead to a speed of the gauge boson proportional to go /ml. If we expect 
this gauge boson to represent a physical photon, then we need this speed 
to be of the order of the speed of light in vacuo (i.e. unity in our system 
of units). This implies that we have to be in a non-perturbative regime, 
where go is of order of ml, and, thus, the one- loop analysis is no longer vahd. 
For this reason, a realistic study of the IR regime of this theory necessitates 
going beyond the one-loop approximation, which lies outside the scope of the 
current article, and thus is left for a future work. 

5 Conclusions and Outlook 

In this work we have examined the infrared limit of a Lifshitz-type Lorentz- 
violating Yukawa model, in (3+l)-dimensional Minkowski space time, where 
only higher order derivatives are present in the model. This differentiates 
the model, from existing ones in the literature, e.g. [9], where the Lorentz- 
violating terms co-exist, as small corrections, with the standard Lorentz- 
invariant ones. In our case, it is demonstrated that quantum corrections 
restore relativistic-like kinematics in the low energy-momentum regime (IR) 
of the model. In this way one can talk about a dynamically emergent Lorentz 
Invariance in the model, where each particle sees a different effective light 
cone. The group velocities of the fermions are found much smaller than those 
of the scalar field, although both velocities are much smaller compared to the 
speed of light in vacuo, consistent with the IR limit, where the analysis is 
performed. 

Moreover, we have discussed a detailed mechanism for dynamical mass 
generation in the model in the same regime of low-energies. It is demon- 
strated that there is no massless limit in this model, in the sense that the 
quantum-fiuctuations dressed masses never vanish, even in the limit where 




9o .2 



(55) 
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the bare masses go to zero. The model exhibits infrared stabihty, in the 
sense that the zero-bare-mass hmit corresponds to infrared-divergence-free 
resummed quantum corrections. In our approach the resummation is pro- 
vided by a differential Schwinger-Dyson type of approach, in which the rele- 
vant dressed parameters of the model appear in a set of coupled differential 
functional equations. We solved numerically these equations for the dressed 
mass and Yukawa coupling, and arrived at the aforementioned results on the 
IR stability and non-vanishing dynamically generated masses. 

We also discussed briefly an extension of the model, which involves cou- 
pling of the fermions to an Abelian gauge field. The induced dispersion 
relations for the gauge field, which remains massless if we want to preserve 
the gauge invariance, are such that its velocity is also found much smaller 
than the speed of light in vacuo, thereby hampering any interpretation of 
such gauge fields as the physical photons. 

In view of the above properties, a natural question arises as to what 
physical systems, if any, such classes of models might correspond to. In view 
of the emergent Lorentz symmetry, one might view our model as a toy model 
for explaining a rather microscopic origin of Lorentz symmetry, in the low- 
energy limit, in analogy with the case of condensed matter systems, with 
nodes in their Fermi surface, where linearization about them leads to low- 
energy relativistic excitations. Such a Lifshitz-type scenario, with emergent 
Lorentz symmetry, could survive the full inclusion of the Standard Model 
group, since the presence of several coupling constants would allow a fine 
tuning of parameters of the model, in such a way that the effective light cone 
seen by different species of particles could be made identical. 

Another possible application of the (gauged) Yukawa model studied here, 
would be its association with the continuum limit of some Lattice mod- 
els of potential relevance to some condensed matter systems in four space- 
time dimensions, in analogy with the case of planar high temperature super- 
conductivity, where the gauge fields represent fractional statistics in (2-1-1)- 
dimensions. In such models, the gauge field is also different from the real 
photon, as in the gauged model of section 4. Thus, one may view the neutral 
scalar as a (spinless) phonon excitation, e.g. due to lattice ion vibrations, 
while the charged Dirac fermion could represent a spin mode coupled to an 
appropriate gauge potential, which, however, is not the physical photon. In 
condensed matter systems, the relativistic fermions may be associated with 
excitations near a node in the Fermi surface of the microscopic model [10]. 
From our analysis in section 4, one may assume that the massless gauge 
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boson group velocity may represent the velocity of the node of the Fermi 
surface about which we linearized the continuum-model. In such the 
physical significance of the gauge potential may be similar to that of a (3+1)- 
dimensional statistics-changing field for the fermion excitations, which could 
represent some kind of hole excitations, due to doping. 
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Appendix A: One loop effective theory 
One-loop scalar self energy 

For external momentum k and vanishing external frequency, the one-loop 
scalar self energy is, after a Wick rotation. 



-igl tr 



X 



/ 



du! d^p iu!^^ — P^(p • 7) + 
2^(27r)3 a;2 + p6 + ^6 

ia;7° - (p + k)2((p + k) • 7) + 



a;2 + (p + k)6 + ml 
The scalar mass correction is given by the zeroth order in k, which is 



(56) 





f duj d^p + 


— mg 


/ 27r(27r)3(a;2+p6. 


f m6)2 








j 


[f + m6)3/2 






/A\ ln2-l" 

In — + -— 

Vmo/ 3 


+ C»(A 



(57) 



and we finally obtain for the scalar mass 



9l 



InfAVi^ 
\mo/ 3 



(58) 
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The term quadratic in spatial derivatives is given by the quadratic-order term 
in the Taylor expansion of the self energy (1561) in powers of k: 

^(2) . 2 fd^ d^P rA2 + 2p2(p.k)2 8(p-k)^ , , 



where V = u"^ + + m^. Using then the following identity, valid for any 
function /, 

"^'^ (P ■ k)V(p^) = 7^ / dp pV(p^), (60) 



we obtain 



^(2-^ wlk'^ f , 2 f , f 16p^ 6pi° + 14m^p^ 24mV° 

= o dp p du { -TTT-r + 



J J V 31)2 p3 p4 

' ' +777^^:^^ • (61) 



TT 



2m2yo (1 + X6)3/2V3 4(1 + X6) ' 2(1 + ^6)^ 



The latter result has to be identified with iX^^^k'^, which leads to the expres- 
sion (fl0|) for X^^\ with a negative sign. 

One-loop fermion self energy 

The one-loop fermion self energy is, for external momentum k and vanishing 
external frequency, after a Wick rotation, 

y _ ■ 2 f d^p du a;7° - p^(p • 7) + "^o 1 .^^s 

^ ~ y (27r)3 27r V ^2 + (p _ k)6 + ^6 • l^^^J 

The correction to the fermion mass is obtained from the zeroth-order term 
in k: 

,2 r ^3 



f (27r)4 7 ^ 1)2 

ig^m^ f p^dp 



7r2 J (p6 + rn6)3/2' 



(63) 
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which yields 

The term hnear in the spatial derivatives is given by the first-order term in 
the Taylor expansion in powers of k: 

^/ - "(2^y^p^'^ — — 



(k-7) 



'''''■(^■7)f7j^,- (65) 



This expression has to be identified with — iA^^^(k • 7), which leads to the 
expression 



2tt^' "J 

167r2 ^ 'V (p6 + m^)5/2- 
(1) 

One-loop vertex 

We show here that the one-loop correction to the coupling vanishes. This 
correction is given by the following three-point graph, for vanishing incoming 
momentum: 

9'" = 9o + 'dt' f^S^, " P^'P • ^ ■ (67) 
" »o 90 y 2x(27r)» (a,2_p8_„6 + fe)3 ^ I 

and, after a Wick rotation, 

^ = go + — duj p dp 

71-^ J [u^ + p^ + rriQy 

= 9o + -r^ P dp' 



The last integral vanishes, and gi^^^ = 510- 
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Appendix B: Evolution equations for the scalar 
potentials 

Taking into account the gradient expansion fl36l) for the constant-field case 
(l)ci 00 and ipciipd ^0 leads to 



5^T 



7^0; - /(p . 7) + [/ 5\p + g) ^ T^'y{p + q) 



^'oo -p\p.^)-U 5\p + q)^ V^h\p + q) 



-U'^j\p + q)^T^p\p + q) 



6(f)p6lPg 

We have to calculate the inverse of the matrix 

/ ri?) T^-^ \ 

r(2) ^ 



(68) 



p(2) p(22 p(2 

V-V" 5/,,/, l/xA 



p(2) p(2_) ^(2) 

\ 'P'^ ipip H 



(69) 



that we find by using the following expansion, valid for any two 3x3 matrices 

{A + i,B)-^ = - A-'BA-% + Oieo). (70) 
In our case, this expansion gives for the elements of the inverse matrix 



1 



^2 _ p6 _ y// 



(71) 



(u;2 - / - y")2 

/l+eo/2 + 0(eo 



X 



2[UfU 



(70a;-p2(p.^))2_f/2_ 



+ (^(e, 
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and 



rr(2)izi = 1 



(72) 



+ 



-2UUfU 



^2_p6_ yn 



(70w-p2(p.^))2_f/2 
= 13 + ^0/4 + 0(^0), 

After a Wick rotation, the integration over (cj, p) gives, for large values of A: 



Tr{/i} 
Tr{/2} 



-iV, / 2A3 
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127r2 \[V"]^ 

iV U" i[U'fU 
2A^W ^ 67r2 

, /2A 



(73) 



In 



U 



[v"]h / {V" - f/2)2 2(y" - m)v" 
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Tr{/4} = V- 



-t\U'fU^ 



In 



iV [U'f , 
+ „}.. In 



[-^"J^ / (V"-f/2)2 2(V"-t/2)[/2 



127r2 (r" - U^) \ [V"- 



1 I , 

2 



where V is the space time volume. The reader should notice that this volume 
factor cancels out in eq. fl35|) . since it appears on both sides of the equation. 



References 

[1] M. Visser, Phys. Rev. D 80 (2009) 025011 [arXiv: 09 02.0590l [hep- 
th]]; B. Chen and Q. G. Huang, Phys. Lett. B 683 (2010) 108 
[arXiv:0904.4565 [hep-th]]; S. R. Das and G. Murthy, Phys. Rev. 
D 80 (2009) 065006 [arXiv:0906.3261 [hep-th]], arXiv:0909.3064 [hep- 
th]; W. Chao, arXiv:0911.4709 [hep-th]; R. lengo and M. Serone, 



25 



[2: 

[3] 

[4] 

[5] 
[6] 

[7 

[8 
[9 



I arXiv: 1003.44301 [hep-th]; G. Koutsoumbas and P. Pasipoularides, 
I arXiv: 1 006.3199 [hep-th]; J. M. Romero, J. A. Santiago, O. Gonzalez- 
Gaxiola and A. Zamora, "arXiv:1006.0956 [hep-th]; K. Anagnostopoulos, 
K. Farakos, P. Pasipoularides and A. Tsapahs, arXiv: 1007.03 55 [hep-th]. 

A. Dhar, G. Mandal and S. R. Wadia, Phys. Rev. D 80 (2009) 105018 
[arXiv:0905.2928 [hep-th]]; A. Dhar, G. Mandal and P. Nag, Phys. Rev. 
D 81 (2010) 085005 [a rXiv:0911.5316 [hep-th]]. 

J. Alexandre, K. Farakos and A. Tsapalis, Phys. Rev. D 81 (2010) 
105029 [ar Xiv: 1004.42 01 [hep-th]]. 

J. Alexandre, K. Farakos, P. Pasipoularides and A. Tsapalis, Phys. Rev. 
D 81 (2010) 045002 [a rXiv:0909.3719 [hep-th]]. 

P. Horava, Phys. Rev. D 79 (2009) 084008 {arXiv:0901. 37751 [hep-th]]. 

R. lengo, J. G. Russo and M. Serone, JHEP 0911 (2009) 020 
|arXiv:0906.3477l [hep-th]]. 

J. Alexandre and J. Polonyi, Annals Phys. 288 (2001) 37 
|arXiv:hep-th/0010128|. 



J. Alexandre, Annals Phys. 312 (2004) 273 [ arXiv:hep-th/0310093| . 

D. Anselmi, Phys. Rev. D 79 (2009) 025017 [ arXiv:0808.3"475] [hep-ph]]: 
Eur. Phys. J. C 65 (2010) 523 [arXiv:0904ji49| [hep-ph]]. 



[10] For a partial list of references see: G. Baskaran and P. W. An- 
derson, Phys. Rev. B 37 (1988) 580; R. B. Laughlin, Science 242 
(1988) 525. N. Dorey and N. E. Mavromatos, Phys. Lett. B 250 
(1990) 107; Nucl. Phys. B 386 (1992) 614; G. W. Semenoff and 
N. Weiss, Phys. Lett. B 250 (1990) 117; A. Kovner and B. Rosen- 
stein, Phys. Rev. B 42 (1990) 4748; J. Polchinski, |arXiv:h ep-th/9210046[ 



K. Farakos and N. E. M avromatos, Phys. Rev. B 57 (1998) 3017 
|arXiv:cond-mat/9611072]; L. Balents, M. P. A. Fisher and C. Nayak, 
[arXivxond- mat/9803086, N. E. Mavromatos and S. Sarkar, Phys. 
Rev. B 62 (1999) 3438 | arXi v: cond-mat /99 1232 3]; J. Alexandre, 
N. E. Mavromatos and S. Sarkar, Int. J. Mod. Phys. B 17 (2003) 



2359 arXiv:cond-mat/0303611|; N. E. Mavromatos and J. Papavassil- 



iou, arXiv:cond-mat/0311421, I. F. Herbut, Phys. Rev. Lett. 94 (2005) 



26 



237001 |arXiv:cond-mat/0410 557]; S. Sachdev, M. A. Metlitski, Y. Qi 
and C. Xu, Phys. Rev. B 80 (2009) 155129. 



27 



